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Abstract. 

We construct a phenomenological model of electroweak-scale inflation that is in 
accordance with recent cosmic microwave background observations by WMAP, while 
setting the stage for a zero-temperature electroweak transition as assumed in recent 
models of baryogenesis. We find that the scalar spectral index especially poses 
tight constraints for low-scale inflation models. The inflaton-Higgs coupling leads 
to substantial mixing of the scalar degrees of freedom. Two types of scalar particles 
emerge with decay widths similar to that of the Standard Model Higgs particle. 
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1. Introduction 

With the results from the WMAP mission it has become relevant to critically review 
models of inflation, especially with regard to the scalar spectral index. While clearly still 
susceptible to improvement, the cosmic microwave background (CMB) observations are 
accurate enough to rule out certain (classes of) models, see e.g. jSHHEJ- In this paper 
we consider electroweak-scale inflation, which turns out to be indeed tightly constrained 
by the spectral index. 

The motivation for looking at electroweak-scale (i.e. of order 100 GeV) inflation 
is twofold. Firstly, it is interesting to see if one can construct a working model of 
inflation with just minimal extensions of the Standard Model (SM) of particle physics, 
and to derive what kind of additional constraints such a coupling to the SM puts on 
an inflation model. The second (main) motivation has to do with baryogenesis, the 
production of the observed baryon asymmetry in the universe. As reviewed in jfjj. all 
necessary ingredients for baryogenesis (baryon number violation, C and CP violation, 
and non-equilibrium) are present in the SM. This provides a strong motivation for trying 
to construct a working model of electroweak baryogenesis. However, the current lore is 
that in a standard finite-temperature electroweak transition both the CP-violating and 
the non-equilibrium effects are too small to be able to account for the observed baryon 
asymmetry. These problems may be resolved in the context of tachyonic preheating at 
the electroweak scale after low-scale inflation [TUBUSES El E] A tachyonic electroweak 
transition is strongly out of equilibrium, and the fact that the process takes place at 
zero temperature at the end of inflation may maximize the effectiveness of CP violation 
|13j . In addition, the low reheating temperature prevents sphaleron wash-out of the 
baryon number produced 

In this context it becomes important to check that the models that combine low- 
scale inflation with tachyonic electroweak preheating satisfy all the new observational 
constraints from WMAP. Low-scale inflation has been considered in many papers (see 
for instance [HI HH1 EH EE1 HZ] ) • In this paper we build in particular on [10 , which was 
also motivated by the problem of electroweak baryogenesis. The main idea is that we 
have a kind of hybrid inflation model ^Hl in which inflation is driven by a nearly 
constant potential energy of order (100 GeV) 4 [7JIS1, while one field (the inflaton) slowly 
rolls down its potential and the other field (the Higgs) is in a local minimum at zero. 
Once the inflaton passes a critical value, the local minimum for the Higgs field develops 
into a local maximum and both fields roll down rapidly to the absolute minimum at 
a non-zero value of the Higgs field, thus breaking the electroweak symmetry. As was 
shown in jTHj, ordinary hybrid inflation models in which the inflaton rolls from large 
field values towards zero are not viable at low scales because of large quantum loop 
corrections. This problem can be avoided in inverted hybrid inflation models, in which 
the inflaton rolls away from zero and inflation takes place at very small field values. 
Note that unlike standard hybrid inflation, slow-roll inflation ends in this case before 
the critical value is reached, instead of the end being caused by the phase transition, 
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so that the slow-roll inflation stage and the phase transition can be considered as two 
separate processes. 

The paper (THj was written before WMAP, and the authors did not study the 
spectral index. As we will show in this paper, their model gives a spectral index that is 
too low according to WMAP. In somewhat more general low-scale inflation models 
were considered, although not from the point of view of electroweak baryogenesis, but 
these models still appear to be incompatible with WMAP. We shall show that one can 
improve these models to obtain a spectral index that lies comfortably within the WMAP 
confidence levels. 

At first sight it may seem that one can always fine-tune a model with sufficient 
parameters to satisfy the constraints, but this is not necessarily the case within a set 
of reasonable rules. To formulate these rules we start from the point of view that we 
are constructing a purely phenomenological model, a minimal extension of the SM that 
introduces only one extra inflaton field in order to describe the history of the universe 
during and after inflation. We stress that there is nothing wrong with fine-tuned 
parameters in a phenomenological model, as the phenomenologically very successful 
SM shows. In incorporating a slow-roll inflationary regime compatible with the CMB 
measurements, one is led to an inflaton potential with non-renormalizable couplings. To 
constrain this potential we assume a polynomial approximation, such that there is only 
a limited number of terms to be parametrized. 

This leads to a tight fit when we also incorporate the scenario of tachyonic 
electroweak baryogenesis, which requires that the inflaton field ends up far from the 
slow-roll regime with a vacuum expectation value similar to that of the Higgs field. The 
inflaton is a gauge singlet and it couples only to the radial (gauge- invariant) mode of 
the Higgs field. This coupling should induce a sufficiently fast tachyonic electroweak 
transition to make baryogenesis possible, without being unrealistically large. It implies 
a considerable mixing between the inflaton and Higgs modes, and the model predicts 
the existence of (only) two scalar particle species with electroweak-scale masses. Up 
to mixing-angle factors, their decay widths are similar to that of a SM Higgs with the 
same mass. The model should therefore be falsifiable by accelerator experiments, in 
particular with the LHC. 

An important issue with any slow-roll inflation model is the question whether 
the assumed flatness of the effective potential is consistent with basic properties of 
quantum fields, with 'quantum corrections'. We investigate this by calculating one-loop 
corrections to the effective potential. This exercise also led us to a rough estimate 
of the scale at which the model may be expected to break down because of its non- 
renormalizable and strong couplings. 

The outline of the paper is as follows. In section |21 we first address the number of 
e-folds of inflation between horizon crossing of a WMAP-observable scale and the end 
of inflation, which number is crucial for the computation of CMB observables. Contrary 
to the generic situation, there is little uncertainty here because the (p) reheating of the 
universe and the onset of the radiation-dominated era are reasonably well understood 
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in this model. Next, in sectional we review the model of and show that its spectral 
index disagrees with WMAP. The implied inflaton-Higgs mixing is studied in section 0] 
We then show in the following section (plus appendix A) that, by adding two additional 
terms to the potential and tuning the coupling parameters to a certain extent, values 
for the scalar spectral index in agreement with WMAP can be obtained. In section |H1 
(plus appendix B) we calculate one-loop quantum corrections to the effective potential 
and study the implications. Finally, section [7| summarizes our conclusions. 



2. Number of e- folds 



One of the most important differences between low-scale inflation and 'normal' inflation 
taking place around the GUT scale is that the number of e-folds of inflation between 
horizon crossing of the observationally relevant modes k and the end of inflation is much 
lower. An expression for Nk is derived as follows [20] (see also the recent papers [2"Hl2"2"]): 

k dH-f^H Or a e a re h a eq i^H _AT, ( PcqV^/i . \_1 Kp\i /i\ 

- - - ' ' * (1 + z cq ) . (1) 



^reh ^eq VPreh/ 1 VSH 

Here a is the scale factor, H = a/a the Hubble rate, z is the redshift, the subscripts 
0, H, e, reh and eq denote evaluation now, at horizon crossing (k = aH), at the end 
of inflation, at the end of reheating and at radiation-matter equality, respectively, and 
k is the inverse reduced Planck mass, k 2 = 8ttG (k -1 = 2.436 x 10 18 GeV). Here we 
used the fact that p oc a~ 4 during radiation domination and the Friedmann equation 
to rewrite Hr. Furthermore we made use of the fact that electroweak (p)reheating is 
nearly instantaneous on the Hubble time scale at the end of inflation, a e /a re h = 1, since 
its time scale is of order of 1 GeV -1 for the SM degrees of freedom j2S], whereas the 
Hubble time H~ x at the end of electroweak-scale inflation is of order 10 14 GeV -1 . At 
radiation-matter equality the energy density p eq is twice that in non-relativistic matter, 
p eq = 2Q m (l + z cq ) 3 p c0 , with p c0 the critical density at present, p c0 = 3k~ 2 Hq. Moreover, 
in our model p re h, p e and pu are all practically equal, and so we find 

where we used P fi m = 0.29, 1 + z COi = 3455, H = 73 km/s Mpc -1 , the WMAP pivot 
scale of k/ao = 0.05 Mpc -1 and an inflationary energy scale of p\[ A = 100 GeV. Hence, 
this number of e-folds Nk is much smaller than the 50-60 one gets in the customary 
models where inflation takes place at much higher energy scales. 

As we will show below, the scalar spectral index n = n— 1 is approximately inversely 
proportional to iV^. This means that the smaller N^ of low-scale inflation makes it more 
difficult to satisfy the WMAP constraint that n should be close to zero. More precisely 
the constraints from WMAP (including CBI and ACBAR, but no other experiments) 
for the amplitude |5k| 2 and spectral index h of the CMB power spectrum are given by 
PP (converted to our normalization, see the definitions in (JHJ)): 

|5 k | 2 = (3.8 ±0.5) x 10 -10 , n = -0.03 ±0.03. (3) 
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3. Core model 

3.1. Model and WMAP constraints 

The model proposed in contains a scalar field a (the inflaton) in addition to the SM 
Higgs field 0. The effective potential of the scalar fields has the formf 

V(a, <P) = V Q - ^a p a p + -^a q a q - \\ a4> o 2 tf + ^/i 2 2 + ^A (0 2 ) 2 , (4) 



with integer q > p > 2. The authors of [10] arrived at the values p = 5, q = 6. As we 
shall explain below, the value of p is fixed by matching the inflationary (small a) part of 
the potential to the SM physics part where a and are near their vacuum expectation 
values. First we need to establish the connection with the CMB data (J5J). 

We choose the inflationary energy scale = 100 GeV. This choice guarantees 
that after preheating and thermalization, the temperature T re h is substantially below 
the electroweak crossover temperature T c ~ 70 GeV [21], thereby avoiding sphaleron 
washout of the generated baryon asymmetry. The reheating temperature can be 
estimated as T re h = [30V A o/(vr 2 ^)] 1 / 4 ~ 0.43 V 1/4 , with g* = 86.25 the effective number 
of SM degrees of freedom below the W mass.§ 

Initial conditions are assumed such that the inflaton has a tiny but non-zero value 
do ~ 10~ 10 GeV. The Higgs field is assumed to be in the ground state corresponding 
to this value of a, i.e. = 0. Only when a reaches the critical value o~ c = /V a/A<t</>, 
which happens long after inflation has ended in this model, will roll away from zero 
and break the electroweak symmetry. This means that we can consider the single-field 
slow-roll inflation stage and the phase transition as two separate processes. 

At the tiny values of a relevant for inflation the a q term in the potential is irrelevant. 
Taking only the first two terms into account, the model can be solved analytically in 
the slow-roll approximation in terms of the number of e-folds N since the beginning of 
inflation: 

a ,N — ~7TEF) = ~ TTFT - <t{N)= [<t p + {2-p)-^-N 



3H 2 k 2 V y ' \ u k 2 V 

2 

1 2 / \2 a p 2p-2 ~ a ,NN ~ , 1 x OL p 2 , . 

e =o K ( a ,N) = 7H772 a > V = ~ e={p-l)^f-a p , 5 

where e = —H/H 2 and fj = a /(Ha) are the slow-roll functions, and we have used 
the fact that e <C fj, because a p a p V during inflation, to neglect the e term in the 
expression for fj, as well as to set H 2 = (k 2 /3)Vq. Defining fj = 1 as the end of inflation, 
we can now compute the scalar amplitude and spectral index to leading order in slow 
roll (see e.g. [201123): 

\ Here <f) 2 /2 stands for (p^p, with <p the usual complex SU(2) Higgs doublet of the SM. 
§ If there are in addition three thcrmalizcd relativistic sterile neutrinos, = 91.5. 
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(GeV) 
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V = 


4 


2.71 x 1(T 12 


1.4 x 1(T 9 


3.6 x 1(T 10 


1.9 x 10" 


-56 


p = 


5 


1.51 x 1(T 3 GeV" 1 


1.4 x 10~ 9 


5.4 x lO" 10 


1.0 x 10" 


-56 


p = 


6 


6.75 x 10+ 5 GeV" 2 


1.5 x 1(T 9 


7.2 x 1(T 10 


7.5 x 10" 


-57 



Table 1. Some quantities in the inflation model with potential V = Vq — ^a p a p 
with Vq = (100 GeV) for p = 4,5,6, as follow from applying the WMAP amplitude 
constraint for the mode that crossed the horizon at Nn, defined as Nk = 22 e- folds 
before the end of inflation N e . 

fi = -2% - 4?h = -2 - « _l_Pjzl (6) 

' (p-l) + (p-2)N k Nkp-2' y) 

using again that en <C fju- Hence we see that h is indeed approximately inversely 
proportional to N k . 

From © we see that the best (least negative) value one can get is n = —0.084 (in 
the limit of large p). Actually the situation is even worse than this, for two reasons. 
In the first place we find in numerical studies that generically there are about 2 more 
e-folds of inflation after fj = 1 has been reached, which means that one should use 
Nk ~ 22 instead of 23.8. Secondly, one cannot take an arbitrarily large p, because, as 
we will show below, p > 5 is not compatible with the constraints from the Higgs sector. 
This means that the upper limit is actually h = —0.114 in this model, which is 3a away 
from the WMAP result (J3J). 

The coefficient a p is determined by fitting the amplitude (JOJ) to the WMAP value 
(JHJ). Results are given in tabled for the cases of p = 4, 5, 6. From this we can derive at 
which value o~ c of the inflaton slow-roll inflation ends (fj — 1), and at which value oh 
horizon crossing of the scale under consideration occurs, see the table. Note that e at 
the end of inflation, also given in the table, is still tiny. 

3.2. Higgs sector 

Next we look at the matching to the Higgs sector, i.e. we consider the stage after inflation 
when a has grown larger than a c and <fi is no longer zero. Here we have the following 
constraints. Denoting the value of the fields in the absolute minimum by v a and v^, 
respectively, we have firstly the two conditions that the first derivatives of the potential 
with respect to a and in the point {v a ,v^) vanish. Secondly, we demand that the 
second derivative with respect to in this point is equal to the the diagonal Higgs mass 
squared, m 2 . The actual particle masses are to be obtained from a diagonalization of 
the cr-0 mass matrix, which we will discuss in section 0] Thirdly, there is the condition 
that l / (f (T ,f 9 i) = 0, so there is no residual vacuum energy (or it is at least negligible 
compared to the electroweak scale). 

Taking v^, and \x as additional input parameters, these four constraints lead to 
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expressions for A CT 0, A^, a q and v a : 

Kef, = (p? + -mjj v~ 2 , A^ = a q = ^/i 2 + ^m 2 ^ v\v- q + a p v p ~ q . (7) 

The vacuum expectation value of the inflaton, v a , cannot be given analytically, but can 
easily be computed numerically from the condition that V(v a , v$) = 0, given the above 
relations. We take 

v^, = 246 GeV, m = 200 GeV, ll = 100 GeV, (8) 

where is the usual SM value. The value for the Higgs mass is of course not known 
yet, but we require that the eigenvalues of the cr-0 mass matrix are above the current 
lower bound on the Higgs particle mass of 114 GeV It turns out that to satisfy 
this lower bound we need to choose a much larger value for the diagonal Higgs mass; 
200 GeV will be sufficient. 

A further condition comes from baryogenesis: the inflaton-Higgs coupling should 
be large enough that the tachyonic transition is sufficiently out of equilibrium. The rate 
of change of the effective Higgs mass squared when a crosses a c = iij v/A^ is determined 
by a dimensionless velocity parameter u: 



1 d{n 2 - A^cx 2 ) _ y/A^ da 
U -<ly? dt ~ y? dV t_tc ' [J) 

where t c is the time when a = a c . We have used /i to set the scale. The value of 
this mass is not critical here (although it may influence the precise value of the baryon 
asymmetry generated) and we have chosen li = 100 GeV. We shall assume that u £ 0.15 
for sufficient baryogenesis. || The value of u can be estimated from energy conservation 

= V o-V(<Tc,0)=cV o , 0<c<l \ a4> = ^—u 2 . (10) 

Requiring u ^ 0.15 means that A^ ^ 0.01/c > 0.01 (for our choices \i = 100 GeV and 
Vq^ = 100 GeV). On the other hand, A^ should not be too large in order to have 
radiative corrections under control, say < 1. 

From the expressions given in ((7|) and (jH)) we find the results given in table El Now 
we can draw the conclusion alluded to before: p = 5 is the only value that satisfies all 
constraints. For p < 4 the coupling A CT ^ between the inflaton and the Higgs is too small 
for baryogenesis, while this case is also worse for the scalar spectral index n. For p > 6, 
A CT9 i is much too large, making quantum corrections uncontrollable. Note that even the 
extreme p = 6 case with q = 7, fi = 1 GeV and = 132 GeV (which means that the 
smallest eigenvalue of the a~4> mass matrix is 115 GeV) gives a that is still much too 
large (X a( p = 625), so that this conclusion does not depend on our parameter choices. 
Hence, we really cannot go beyond p = 5 when we try to maximize the spectral index. 

The basic reason for the limit on p is the fact that p > 5 implies a rapid turning 
down of the potential as a increases and, together with the low value V^ 4 = 100 GeV, 

|| Our u is equivalent to the velocity V of |27l 111) . 
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3.2 x 10" 22 
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1.45 x 10 5 


V = 


■5, 


q = 


6 


0.36 


0.33 


8.4 x 10~ 6 


GeV" 2 


288 


V = 


6. 


q = 


7 


2.0 x 10 3 


0.33 


3.0 x 10+ 5 


GeV" 3 


3.90 


V = 


6, 


q = 


8 


3.1 x 10 3 


0.33 


2.9 x 10+ 5 


GeV" 4 


3.09 



Table 2. Values for the parameters and inflaton vacuum expectation value v a in the 
full model f° r different choices of (p, q) , following from constraints on the Higgs 
sector of the model, as explained in the main text, using the values given in table [U 
and equation 



this leads to a very small v a ; this in turn requires a very large X a( f ) in order to obtain 
the required value A^f 2 — /i 2 = m 2 /2. 

We end this section by noting that for the p = 5, q = 6 case, n = —0.114; 
furthermore o c = 166 GeV, which satisfies a e <^ cr c < v a as we assumed, and c = 0.09 
(cf. dm))). Note that A CT9 j = 0.36 satisfies the requirement X a ^ > 0.01/c. The energies 
corresponding to the values of and a§ are somewhat above the electroweak scale: 
1/ob = 661 GeV, l/y/ai = 344 GeV. 

4. Inflaton— Higgs mixing 

In this section we take a first look at the possibility of testing the model with accelerator 
experiments. The shape of the potential near its absolute minimum determines the 
particle masses and interactions that we can measure in the laboratory. The inflaton 
and Higgs fields have the same quantum numbers after electroweak symmetry breaking, 
and consequently the emerging particles correspond to a mixture of the two fields. The 
particle masses are given by the eigenvalues of the a— (f) mass matrix, 

f \r° 1/^ ) ' ° = v a , <p = (11) 

and for the core model with p = 5, q = 6 they are given by mi = 385 GeV, m2 = 125 
GeV, safely above the current experimental lower bound for the Higgs mass of 114 GeV. 
The mixing angle defined by 

0i = acos£ — 0sin£, 2 = crsin£ + 0cos£, (12) 

where <pi and 02 are the mass eigenmodes, is given by sin£ = 0.43, cos£ = 0.90, or 
f = 0.44. 

Because mi > 2m2, the heavier particle can decay into two lighter ones and it is of 
interest to calculate the decay rate (see e.g. [28J): 
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with 7 12 2 the three-point coupling: 

d 3 V 

We find 7122 = 53 GeV and a decay rate I\^2+2 = 56 MeV. However, the mixing 
into the Higgs proportional to sin £ facilitates decays into the modes allowed for a SM 
Higgs particle with mass mi, e.g. the decay into two W bosons, with a much larger rate 
reaching 100 GeV for a mass ~ 500 GeV [26J. Since the latter decay mode is forbidden 
for the lower mass particle 2, it decays predominantly via the bb channel with a much 
smaller rate of order 10 MeV j!26]. Hence the branching ratio via 1 — > 2 + 2 will be very 
small. So the model predicts 

ri_oc ^ sin 2 £ r H ^ x , m (j) = m 1 , (15) 
T 2 ^x ^ cos 2 £ r H ^ x , "V = m 2 , (16) 

for a dominant decay mode X of the SM Higgs particle. 

Let us conclude this section by discussing the number of parameters of the 
model. Consider coupling the gauge-singlet inflaton field to the SM Higgs field in a 
renormalizable fashion. To start out with this would introduce three new parameters 
beyond those of the SM: 

ml = V aa (v a , v<f>), AjE - V aaaa (v a , Vj) (17) 

and X a( p. Lacking the symmetry a — ► — a, we should also include V^^^iy^.v^) and 
V,o<f><t>{ v o-i v 4>)- The linear term in the expansion around the minimum vanishes since 
V a (v a , v$) = 0, which determines v a . We should also keep in mind the cosmological 
constant, which is needed to control the energy density of the vacuum (approximated 
by zero in this paper) and which is usually not included in the parameter set of the 
SM. So within the renormalizable class of models we should not be surprised to find 
seven new parameters. The scenario for baryogenesis suggests the introduction of two 
more parameters in order to be able to adjust the height and the first derivative of 
the potential at the spinodal point a c \ V(cr c ,0) and V a (a c ,0) (where a c is determined 
by V^(a c ,0) = 0). These two conditions suggest that we need more freedom in the 
potential, e.g. provided by the coefficients of non-renormalizable a 5 and a 6 terms. 

In the core model all these parameters are basically set by the inflationary 
parameters Vq and a p , except for the cosmological constant which is controlled by a q . 
It would of course be surprising if the parameters thus obtained were just right for the 
phenomenology in the electroweak domain. To influence their values (e.g. to control 
ml and hence also m\) we would have to add further terms to the potential in the 
electroweak regime. For the case p = 5, q = 6 these would be a 7 , a 8 , .... Experiment 
will inform us of the need of such terms. Moreover, in the inflationary domain there 
are terms that have been arbitrarily set to zero in the potential of the core model with 
p — 5, namely the a 2 , a 3 and a 4 terms. These will be considered in the next section. 
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n-l 



-0.08 




Figure 1. The scalar spectral index h = n— 1 as a function of the dimensionless mass 
parameter x = a2/(« 2 Vo) in the model (|18fl . as given in equation H20|) . with Nf. = 22 
and plotted for p — 4 (solid), p = 5 (long dashes) and p = 6 (short dashes). 



5. Improving the core model 



5.1. Adding a quadratic term 

To bring the scalar spectral index closer to the central WMAP value we add a negative 
quadratic mass term for the inflaton to the potential, as in the class of models studied 
in 



V(a,(j)) = Kore( 



1 2 



with V core (a,<j)) defined in (J3J). In this case the slow-roll system (considering only 
Vq — \ct2O~ 2 — -a p a p ) can still be solved analytically: 



a(N) 



a 2 



2-p 



Op 

a 2 



On 



2k 2 V£ 



a 



1 + 



a. 



P a P-2 



exp (2-p) 



7/ 



OL2 
Q'2 



N 



1 

2-p 



(p-1 



(19) 



a 2 J n 2 V \ a 2 

which agrees with (j^J) in the limit that a 2 — > 0. 

In a derivation analogous to the one in the previous section, we find for the spectral 
index the following expression: 

(p-2)(l -x) + (1 + (p-2)x)e (p - 2)xNk .„ o 2 



-2a; 



with a; = 



(20) 



(ar - 1) + (1 + (p - 2)a;)e(p- 2 ) :r ^ K 2 y 
The curve n(x) is plotted in figure ^ for p = 4, 5, 6 and Nk = 22. We see that the 
situation can be somewhat improved compared to the massless (x = 0) case by choosing 
a 2 such that x ~ 0.023, i.e. a/q^ = 6.3 x 10~ 16 GeV, which for p = 5 changes n from 
—0.114 to —0.092. However, this is still a 2a deviation from the central WMAP result. 



Note that the added quadratic term is completely negligible for a 



and its only 
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effect for the Higgs sector comes from changing the value of a p that follows from the 
WMAP amplitude (to a$ = 5.0 x 10 -3 GeV -1 for the case p — 5). This does not change 
the conclusion about p > 5 being excluded. 

5.2. ... and a quartic term 

Including in addition a quartic term for the inflaton allows for more variation in the 
spectral index: 

V(a, 4>) = V COTe (a, 4>) - ^a 2 cr 2 + ^a^cr 4 , (21) 

with V core (a, (ft) defined in (jlj), where we have now fixed p = 5, q = 6. This means that 
for the inflationary part we consider the potential 

Mnfl(^) = V - \a 2 (J 2 + - la 5 a 5 . (22) 

2 4 5 

For the analysis it turns out to be useful to make the following definitions: 
TJa 2 a\ a 2 _ Sa 5 

R = —r^ X = ^v' a = ^r a (23) 

(this is the same definition for x as in the previous subsection). Assuming that all 
a's are positive there is always a maximum at a = 0, and for a > there are three 
different cases depending on the value of R. For R > 27/4 the potential has a negative 
second derivative for all positive a, for 4 < R < 27/4 the second derivative changes sign 
from negative to positive and back to negative, and for R < 4 there is an additional 
minimum-maximum pair. For R = 4 the potential has a flat plateau around a = 2. 
We restrict ourselves to the region with R > 4, because an additional minimum would 
stop the inflaton from ever reaching the critical value a c (ignoring quantum tunnelling 
which we do not consider in this paper) and the electroweak symmetry would not be 
broken (nor would inflation end at all). In the full model there is of course no drop off 
to minus infinity, because after inflation the \a^a & term will come into play to create 
the absolute minimum at v a . 

Using the slow-roll equation of motion, we can derive expressions for e(a) and ry(er) 
in this model as well, leading to the following results for the amplitude and spectral 
index, still in terms of the inflaton field: 

N 2 = ^2 (f ) 3 - 34 + 4)" 2 , n = -| (R - 94 + 44) • (24) 

To determine <7h we need to solve the equation of motion. Although it can be integrated 
analytically to give N(a), the inversion to obtain a(N) has to be carried out numerically. 
The horizon-crossing value <th is then the field value 22 e-folds before the end of inflation, 
which is defined by the relation fj — 1. The analytical results used in this procedure are 
given in appendix A. 

An important result that can be seen from the equations in the appendix, is that 
5h depends on the parameters R and x only. That means that the same is true for n. A 
contour plot of n as a function of x and (R — 4) _1 is given in figure El (we use (R — 4)~ x 
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Figure 2. Contourplot for n = n — 1 as a function of the parameters x and (R — 4) _1 
(defined in l |23| l) in the model l|22|) for the mode that left the horizon Nk = 22 
e-folds before the end of inflation. The contours are (from outside to inside) for 
n = -0.12,-0.09,-0.06, -0.03,0.00, +0.03, +0.06, respectively. This means that the 
red region (between the contours —0.06 and 0.00) corresponds with the la WMAP 
constraint @, while the yellow/green regions (near the axes) indicate values that are 
too low, and the blue regions values that are too high. The blue regions also correspond 
with a blue (n > 0) spectrum. 

for the ordinate to avoid the curves getting squashed in the region R J. 4). We conclude 
that there is a parameter region (the red area) where a spectral index compatible with 
the la WMAP constraint © is produced. For R < 27/4 (i.e. (R - 4)" 1 > 4/11) it is 
even possible to get n > (the blue region). On the other hand, the amplitude does 
not just depend on R and x, but also on a 4 explicitly, as can be seen from (J2U). 

Now we can apply all the constraints and determine the parameters in the following 
way. First we determine x and R from the spectral index constraint. Of course this 
cannot be done uniquely; basically we are free to choose e.g. x (within certain limits as 
indicated in figure EJ) and then R is fixed by the constraint. Actually there are still two 
possibilities for R, but it turns out that the smaller value (corresponding to the upper 
branch in figure EJ) leads to smaller non-renormalizable couplings a 5 and a^, and also 
to a smaller dimensionless coupling A CT (cf. (jl7|0 . which we consider more acceptable. 
Next «4 is fixed by the CMB amplitude constraint. Finally, the other parameters are 
determined by the coupling to the Higgs sector, in a way completely analogous to the 
treatment in sections I3.2l andl4l Setting n = —0.03, the central WMAP value, an explicit 
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Figure 3. Variation of the dimensionful (top) and dimcnsionless (bottom) parameters 
in the model described in section T5. 21 as a function of the one free parameter x when 
all CMB (with h = —0.03) and SM constraints are applied (and choosing the smaller 
value for R). 



example with x chosen to be 0.1 is: 

£ = 0.10: R = 4.32, ouj 2 = 1-30 x 10~ 15 GeV, a 4 = 4.36 x 10~ 12 , 
a- 1 = 357 GeV, a~ 1/2 = 238 GeV, = 0.46, 

A = 0.33, X a = 8.6, = 246 GeV, v a = 255 GeV, 
mi = 427 GeV, m 2 = 128 GeV, f = 0.39, 7 i22 = 65.5 GeV. (25) 

For a more general picture see figure El where all the parameters have been plotted as 
a function of x (and choosing the smaller value for R). We see that as x increases, the 
non-renormalizable mass scales a^ 1 and ctg also increase and that the (rather large) 
inflaton self-coupling X a decreases. Since smaller couplings are generically easier to deal 
with, the larger x values would be favoured, but we shall not pursue this aspect further. 
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6. Loop corrections 

Thus far we have assumed the potential V to be an effective potential that approximately 
describes some extension of the SM, i.e. including all quantum effects. Nevertheless, 
since the calculation of the spectral index is based on the growth of quantum fluctuations 
during inflation, it is important to try to ascertain that the back-reaction of these 
fluctuations on the inflaton and on the SM is under control. After all, the a values vary 
over some twelve orders of magnitude from the inflationary to the electroweak domain. 
We therefore investigate in this section one-loop corrections based on V. 

Corrections on the effective potential in the inflationary domain due to a Higgs 
loop have been investigated in [TO], but not those coming from an inflaton loop. There 
is reason for concern, because the four-point self-coupling of the inflaton in the 
minimum of the potential, defined in ()17)) . is rather large. For example, X a = 8.6 for 
the case fl23j) . With such large couplings, computing loop 'corrections' is a hazardous 
endeavour and the calculations in this section are aimed at obtaining insight rather that 
getting quantitative results. 

At strong couplings the non-perturbative phenomenon of 'triviality' comes into play: 
as the cut-off used to define the model is raised the renormalized couplings go down, and 
they even vanish in the infinite cut-off limit. Larger couplings imply a smaller maximum 
value of the cut-off, which is then interpreted as a momentum scale where the model 
breaks down (also known as 'the Landau pole'). For a review of the application to the 
SM see [26, 29J. Assuming that this phenomenon applies also here, we shall tentatively 
use it to estimate the maximum cut-off from the large inflaton self-couplings. 

Another reason for expecting such a scale where the model breaks down is the 
fact that it contains non-renormalizable couplings of dimension larger than four. This 
situation may be compared with effective pion models for QCD. These models are 
typically also non-renormalizable (even non-polynomial) and in the first instance only 
valid up to a few hundred MeV.^f Loop corrections may then further extend their 
validity. A well-developed scheme is chiral perturbation theory, in which all perturbative 
infinities are removed by counter-terms, with new physical constants parametrizing the 
corresponding finite parts jSOJ EH] • 

Before continuing in detail we should add the cautionary remark that the usual 
equilibrium effective potential can in principle not be applied blindly to systems out 
of equilibrium; we interpret it as being only indicative of the back-reaction. More 
sophisticated methods (e.g. Hartree or 2PI [321 03] ) are available, but they are a lot 
more complicated. In the inflationary domain the one-loop potential is complex due to 
the tree potential being unstable, and we shall concentrate on its real part. 

Consider the potential 

V (0) = V - \a 2 a 2 + lo 4 a 4 - ±a 5 a 5 + l -a & o* - ^a 2 2 + + \\^\ (26) 

% We recall also the linear sigma model, in which the self-coupling is large, A = m^/(2/^) ~ 20-45 for 
a sigma resonance in the range 600-900 MeV. 



Electron) eak- scale inflation, inflaton-Higgs mixing and the scalar spectral index 



15 



which is to be identified with V at tree level. For simplicity we consider here only one 
real Higgs field to avoid the complications of massless Goldstone bosons, which would 
be absorbed by the W and Z bosons in the full SM. The one-loop contribution to the 
effective potential is given by (see e.g. |29|) 

VM = e(ml) + e(m 2 ) + c.t. (27) 

where e(m 2 ) is the ground-state energy density of a free scalar field with mass m, c.t. 
denotes counter-terms and m\ and m\ are the eigenvalues 

™i,2 = \{ml + ml± ^K-m;) 2 + 4m^) (28) 
of the field-dependent mass matrix 

™ 2 ab = V { °i (a,b)e{a, ( j ) }, (29) 
m 2 = m 2 aa = —a 2 + 3a 4 <7 2 - 4a 5 cx 3 + 5a 6 cr 4 - X a( p(f) 2 , (30) 
m l = m l4> = V 2 ~ K<t>° 2 + 3A 2 , m 2 a(j) = -2X a(j ,a(f). (31) 
The ground-state energy density is given by 

d 3 r> 1 



— ^-V^ + p 2 , (32) 
J\p\<a (2vr) d 2 

where we used a spherical cut-off A on the three-momenta. Alternatively, we can use a 
four-dimensional Euclidean cut-off on the usual log-determinant form lndet(m 2 b + d 2 ), 
which leads to 

e(m 2 )= 1 -[ ^M^+P 2 )- (33) 

2 i|p|<A ( 27r ) 

The first regularization (}3*2*|) has the advantage that time is kept real (and 
continuous), which is conceptually attractive for systems out of equilibrium, whereas 
the regularization (|33J) can be applied only in equilibrium with imaginary and somewhat 
fuzzy time. It turns out that both regularizations give very similar results and for 
definiteness we shall continue with the first one, (|32|). 

Within the spirit of a polynomial parametrization, the counter-terms are supposed 
to be polynomial in a and <fi. It is desirable that they are able to cancel all divergencies 
as A — ► oo: the quartic (oc A 4 ), quadratic (oc {m\ + m^)A 2 ) and logarithmic 
(oc ((mf) 2 + (m2) 2 )lnA) ones. This is possible, since the square root in (|28|) drops 
out of the sum m 2 + m 2 ,, and also (m 2 ) 2 + (m 2 ,) 2 = {ml) 2 + (m 2 ) 2 + 2{m 2 a(j) ) 2 is a 
polynomial in a and 4>. We define 

e f (m 2 ) = dm 2 ) - -L ( 2 A 4 + 2m 2 A 2 - (m 2 ) 2 In ^lA) (34) 
327T J V v J 



2 

{ m 2 ) 2 ln-^- + 0(A- 2 lnA), (35) 

where we introduced the renormalization scale v. Dropping terms with negative powers 
of A, the form (J33|) for e reduces also to the form (J35j) (with a different subtraction). 
We shall continue with the full form (|3*4*j) with e(m 2 ) given in (|32|). without dropping 
negative powers of A. 



64vr 2 v 
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The total potential at one loop is 

v = v m + v a)_ ( 36 ) 

At this stage we have used counter-terms corresponding to all the terms in ()26|) to cancel 
the divergent cut-off dependence, as well as a few more: a 7 , <r 8 , o~ 3 (p 2 and <r 4 2 . The finite 
parts of these counter-terms have been assigned values by the subtraction in (which 
is similar to the minimal subtraction used in perturbative QCD calculations). However, 
we still need at least some of these parameters to impose renormalization conditions 
on the potential. A minimal set of conditions is: (i) the real part of the expansion of 
V in a around a = 0, at = 0, coincides with the original up to and including 
a 5 ; (ii) V vanishes in its absolute minimum at a = v a , <f> = and (Hi) = 246 
GeV, unchanged. Requirement (i) expresses our wish to keep the potential unchanged 
in the inflationary domain, (ii) keeps the cosmological constant zero and (Hi) is needed 
for electroweak phenomenology. Furthermore, in a sensible renormalized perturbation 
scheme it is desirable that the particle masses in loop diagrams are kept at the same 
values as in the tree graph starting point, so we add to our minimal set of conditions: 
(iv) v a , (v) m 2 and (vi) m 2 ^ are unchanged in the minimum of the potential. We can 
impose all these conditions by using only counter-terms according to the parameters 
listed in the tree potential (J2EJ). It turns out that m 2 ^ ends up also very close to its tree 
value, so the particle masses ra\ and m 2 are indeed practically unchanged. The unusual 
combination of two classes of renormalization conditions, (i) in the inflationary region 
and (ii)—(vi) in the 'today's physics' region, makes an interesting problem. 

Summarizing, the above renormalization conditions can be implemented by counter- 
terms of the form 
1 



c.t. 



327T 2 

5 



4A 4 + 2(m\ + m 2 2 )A 2 - ((m\) 2 + (m 2 ) 2 ) In 26 — 



fc =o L ° a > -=<f>=o fc - 

* a " ^ ^ s 20 2 » 4 
+ da 6 — - oAptf,— h d(i — + OA^—, 

which satisfies condition (i) (note that m 2 = m 2 and m\ = m 2 in the inflationary region, 
since = 0). The first thing to check is the coefficient of the a 6 term in the expansion 
around a = 0: it should not be large, as this might spoil the neglect of such a term 
in the inflationary domain. For the Higgs loop this has been done already in jlOj . Its 
contribution to the coefficient of cr 6 /6 is — A^/(327r 2 /i 2 ), which for the case (|23|) is only 
—3.1 x 10 -8 GeV -2 , even (much) smaller (in absolute value) than a 6 = 1.8 x 1CT 5 GeV -2 . 
For the inflaton loop the effects are even smaller because the effective mass m a is tiny 
in the inflationary region; although ln(m 2 /V 2 ) is the logarithm of a very small number, 
the (m 2 ) 2 in front of it makes it negligible. Since the terms with different powers of 
a up to and including a 5 are by construction all approximately equal during inflation, 
a rough estimate of the relative importance of the inflaton loop corrections is given by 
a 2 ln(a 2 /z/ 2 )/(647r 2 a 2( x 2 ) = £>(1(T 12 ). 
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We now impose the renormalization conditions (ii)-(vi) in the minimum of the 
potential. These are five linear equations for the four parameters of the counter-terms 
in the third line of (JBZJ), plus the renormalization-scale parameter ln(z/). Using as an 
example the parameter set (|25|) and choosing A = oo for a start, results in 

5a 6 /a 6 = 0.10, 5\^/\ a4> = -0.011, Sfi 2 /^ = -0.058, SX^/X^ = 0.026, 

z, = 939GeV, m^/m™ = 0.95, A^/Af = 2.07, (38) 

where (0) indicates the tree graph value and (0+1) the value derived from the full 
potential (|36|). The renormalization scale v seems to be somewhat large, but small 
changes in \n{y) are of course amplified in v itself. The finite counter-terms look 
reasonably small. The mixing mass is also close to the tree graph value, in 
accordance with our desire to keep the particle masses in the loop contributions close to 
the tree values. The inflaton self-coupling X a has increased by an unpleasant factor of 
about two, which seems to make the calculation untrustworthy. Another way of treating 
the renormalization to deal with this is discussed in appendix B. 

Since we have removed the divergent cut-off dependence with the counter-terms, 
the renormalized effective potential cannot give us information on the scale where the 
model is expected to break down. For this we look at the bare coupling defined in terms 
of the bare potential + at., by 

A( b ) = (l/6)(y(°) + c.t.) i(T(T(T(T , (39) 

evaluated at its absolute minimum. In a simple 4 model with one real scalar field the 
relation between the bare and renormalized coupling is qualitatively well described by 
the one-loop renormalization-group relation (see e.g. |2"9~1 134j ) 

A (b) = A = A + \ 2 f3t + ■ ■ ■, t = In (A/m) + constant, (40) 
1 — p At 

with (3 = 9/(87r 2 ) the one-loop beta-function coefficient. The bare coupling X^ diverges 
when the denominator vanishes at the position of the 'Landau pole', and numerical 
simulations have indicated that this gives an order of magnitude estimate of the limiting 
A (modulo the uncertainties related to methods of regularization and factors like 2e -1 / 4 
in (|34|): for a review see e.g. |34J). As can be seen from the expansion in (|40|). the 
one-loop value of the bare coupling equals twice the tree value at the pole, and we shall 
use this for an estimate of the maximal cut-off A max : 

A( b VA CT = 2 A = A max . (41) 

Using the results from the renormalization procedure in ()38jl we find that the maximum 
value is already reached near A ~ 600 GeV, for which X a = 12.8 (and v = 702 GeV). 

The loop calculation in this section has given evidence that the flatness of the 
inflationary part of the potential can be kept consistent with quantum corrections. It 
appears that the model has to already break down at a fairly low scale. For the example 
(|2*HJl this could be around 700 GeV or perhaps as low as 300 GeV (see appendix B), but 
for larger x values this scale slowly increases. Above this scale new physical input is 
needed. 
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7. Conclusion 

In this paper we investigated the implications of the WMAP results for low-scale inflation 
and found that there are severe constraints. This is essentially because the proximity 
of the observed scalar spectral index n = n — 1 to zero is hard to reconcile with the 
small number of e-folds between horizon crossing of the observable scales and the end 
of inflation in these models. Working in the context of a phenomenological electroweak- 
scale inflation model that allows for tachyonic electroweak baryogenesis and consists 
of one additional scalar field a coupled to the Standard Model Higgs, we were led to 
further constraints on the inflaton-Higgs potential. However, we found that there is a 
range of parameters compatible with a spectral index close to (and even larger than) 
zero in this model. 

The polynomial approximation together with all the constraints led to the 
conclusion that we need a a 5 term in the potential during inflation (as in JU]). In 
addition a 2 and cr 4 terms are needed, and (or instead of the a 4 ) there might be a a 3 
term, but no powers higher than 5 can be present during inflation. The appearance of 
an odd power (a 5 ) implies a local minimum at small negative values of the inflaton field 
(a = being the value where the potential has a local maximum). A universe with initial 
conditions in this negative region would classically never stop inflating, which is why 
we assumed a small positive initial condition for the inflaton field. However, quantum 
tunnelling might very well make the case of a negative initial condition viable as well. 
The use of a polynomial approximation is quite natural for the inflationary region of the 
potential where the inflaton field is small, but it seems somewhat artificial in the large 
field region where the Higgs field comes into play. The odd and non-renormalizable 
power a 5 in the potential appears to be the price we have to pay for keeping the number 
of parameters limited. It is not excluded of course that the model with its non-symmetric 
potential and non-renormalizable couplings can be embedded satisfactorily in a model 
with more symmetry that is also renormalizable, e.g. a supersymmetric extension of the 
Standard Model. 

In the one-loop calculation we found that the very different renormalization 
conditions in the inflationary and electroweak regimes did not lead to unresolvable 
conflicts. In our results quantum corrections do not disrupt the required flatness of 
the potential in the inflationary region. The non-renormalizable couplings and also 
the relatively strong inflaton self-coupling suggest a breakdown of the model already 
occurring at a fairly low scale, perhaps below 1 TeV, depending on the choice of 
parameters. 

We conclude with the important remark that the phenomenological model arrived 
at here can be falsified experimentally through its conspicuous generic feature: the 
existence of two (and only two) particle species with zero spin and masses around the 
electroweak scale, and couplings to the rest of the SM equal to that of the Higgs up to 
factors related to the mixing angle. 
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Appendix A. Some analytical results for the full model with a quartic term 

Starting from equations ()22j) and (|23|) in section 15^1 we find the following slow-roll field 
equation: 

dN = z 7- (A.l) 

Integrating it from horizon crossing (subscript H) to the end of inflation (subscript e), 
we obtain for R > 4 (with N k = N c - jV H ) : 



N k = — 

6X 



with 



2 + Q + 2Q 2 ! a -l + Q + Q- 1 

m 



l + Q + Q 2 \ y /^-(2 + Q + g-i)or + Q + Q~ 1 + Q 2 + Q~ 2 
\a -l + Q + Q- 1 J Q 3 -l v ' 



(A.2) 



, arctan f -j^ffL ) for a < 1 + \Q + 



and 



arctan ( jff+iq-"^- ) + 7t for a > i + ig + |g- ] 

Q z = ]^[R-2 + ^R^R^j R = 2 + Q 3 + Q- 3 . (A.4) 

The field value at the end of inflation that follows from the condition fj = 1 is for x < 1 
given by 




Appendix B. Another renormalization treatment 

We found in (|38[) that the already large inflaton self-coupling A CT increased by another 
factor of two when applying the renormalization conditions. Of course, de facto, X a is a 
derived coupling: it depends on the other parameters of the model and we may just have 
to accept how it turns out. In principle, further experimental information is needed to 
be able to pin down the non-renormalizable couplings beyond a$ and otQ, similar to what 
is done in chiral perturbation theory. To continue, we can pretend that our tree level 
values for couplings such as A CT are actually phenomenologically correct, i.e. consider 
them as experimental input. Alternatively (but leading to the same treatment) we can 
work from the point of view presented at the end of section 0] that A CT is a primary 
coupling and that the other parameters have to be chosen accordingly. (At this point 
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we have assigned values to the non-renormalizable couplings of the a 7 , a 8 , a 3 4> 2 and 
a 4 (j) 2 terms in a somewhat arbitrary way by the subtraction in (J3*4"j) and the value of the 
renormalization scale v.) So we continue this exploration by imposing that in addition 
to v ff , and the particle masses m 2 , and m 2 ^, also the tree graph value of X a is to 
remain unchanged, and that the other parameters have to be chosen accordingly 

We add to our renormalization conditions: (vii) m 2 ^ = m^? and (viii) X a = X& . 
Since A CT is naturally controlled by the a parameters, we add to (|37jl the finite counter- 
terms 

a 7 a 8 

+ <5a 7 y + <5a 8 y, (B.l) 

still keeping the coefficients of a 3 <p 2 and cr 4 2 set by the value of v. The result of this 
exercise is given by 

5a 6 /a 6 = -0.75, 5X a4> /X a4> = 0.15, 5fi 2 /fx 2 = 0.34, SX^/X^ = 0.092, 

v = 3987 GeV, 5a? 1/3 = 217 GeV, (-5a 8 )" 1/4 = 363 GeV. (B.2) 

We see that the scale of 5a-j and 5as is similar to that of and and the other 
counter-term parameters are larger than before in ()38|) but still reasonable, except that 
v has slipped to a rather large value. We could fix v to, say, 400 GeV and use one of 
the a k (p 2 couplings in its place, but then the resulting counter-term parameters for the 
other couplings are larger. Moreover, we prefer to keep v variable since it adjusts itself 
naturally in such a way that the loop correction vanishes as A —>■ 0. 

When we try to determine the maximal cut-off in this case using the estimate ([41)1 . 
we find that this criterion cannot be applied meaningfully because the bare potential 
quickly becomes unstable (no lower bound) as A is increased from zero. This happens 
already for A between 200 and 300 GeV (where v has come down to ~ 300 GeV). 
Such a low value is also suggested by the fact that 250-350 GeV is the scale of the 
non-renormalizable couplings and a§ in the example ()25|) . 
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